## Order of field elements

from PyM import *


# Irreducible polynomial: x**5-x+1
[F,x] = extension(Zn(5),[1,0,0,0,-1,1])

q = cardinal(F)

r = order(x)

show(r)

show((q-1)//r)


y = sqroot(x,F)

show('y =',y)

show('y**2 =',y**2)

show(order(y))